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\Q ", ABSTRACT 

o ■ 

Aims. We test the hypothesis that SO galaxies are the descendants of fading spirals whose star formation has been shut down, by using the 
£Sj . properties of their globular cluster systems. 

I ■ Methods. We estimate the amount by which the globular cluster specific frequency (number of globular clusters per unit V-band luminosity) is 
£3 enhanced in SOs relative to spirals. If the transformation hypothesis is correct, and no clusters are created or destroyed in the process, then this 
"™* . difference provides a measure of the degree to which the SO's K-band luminosity has faded relative to that of its spiral progenitor, which we can 
OO ' compare with the independent values estimated from stellar population synthesis and the SO Tully-Fisher relation. We also explore whether the 
O^l degree to which the globular cluster specific frequency is enhanced in SOs correlates with the colour of the stellar population, as also predicted 
by this hypothesis in which galaxies become redder as they fade. 

Results. We find that, on average, the globular cluster specific frequency is a factor ~ 3 larger for SOs than for spirals, which can be interpreted 
as meaning that passively-evolving SOs have faded on average by about a factor of three from their spiral progenitors. This value fits remarkably 
well with the predictions of stellar population synthesis calculations, and the offset between the SO and spiral Tully-Fisher relations, where the 
SO V-band relation lies ~ 1.2 magnitudes, or a factor of three, below the spiral relation. We also find that the global colours of SO galaxies are 
strongly correlated with their globular cluster specific frequencies: the redder the stellar population of an SO, the larger its specific frequency, as 
' we might expect if we are catching different SOs at different stages of passively fading and reddening. Comparison to the predictions of stellar 
\^ population synthesis models show that this explanation works quantitatively as well as qualitatively. 

f*"*) . Conclusions. These tests strongly support the hypothesis that SO galaxies were once normal spirals, whose star formation was cut off, pre- 
sumably due to a change of environment. We are now in a position to start to make quantitative measurements of when this life-changing event 
'"^jJ , occurred in different galaxies. 

O Key words, galaxies: formation - galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: spiral - galaxies: star clusters 
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1 . Introduction of halo gas c ould lead to a slower "strangulation" of star for 



or nalo gas could lead to a 
mation jLarson et alJl98oT) . 



| Although this scenario is quite plausible, we need 
■ SO, or lenticular, galaxies live at the crossroads between ellip- some more direct evidence that it has actually occurred, 
tical and spiral galaxies in the traditional Hubble tuning-fork Circumstantial evidence for such a transformation over cos- 
diagram, suggesting that they should play a key role in under- mic timescales comes from observations which indicate that 
standing the morphologies of galaxies. In fact, in some cluster the proportion of SO galaxies is substantially smaller in distant 
environment s, SOs are the single most common type of lumi- clusters than in nearby o nes, while spirals sh ow the opposite 
nous galaxy ( Dressle Jl980l) . so clearly understanding how thev trend (see, for example, Dress ier et alJll997i) . but any infer- 
form and evolve is essential if we wish to have a complete pic- ences drawn from such studies always face the criticism that 
ture of how galaxy m orphology is rel ated to galaxy formation one may not be equating comparable systems at different red- 
and the environment. Dressier] Jl980l) also showed that as the shifts. Further evidence that SOs are fading spirals comes from 
fraction of SOs increases in the densest cluster environments, so studies of the Tully-Fisher relation between luminosity and ro- 
the fraction of spirals decreases, naturally suggesting that SOs tation speed for SO galaxies, which is found to be rather broad 
are simply "dead" spiral galaxies that have had their star for- and offset from the relation for spirals in the sense that the SOs 
mation shut off by their surroundings, and are now quietly fad- are systematically fai nter than the spirals. A recent analysis by 
ing away as their stellar populations age. Indeed, a number of iBedregal et alJ (J2006) found that the SO relation in the fi -band 
mechanisms have been suggested that might cause such a trans- is offset by ~ 1 .4 magnitudes from the ISakai et aD J2000I) re- 
formation: ram-pressure stripping of disk gas could remove th e lation for spirals, with a scatter of ~ 1 magnitude. Such an 
raw material of the next stellar generation ddunn & Gotl 19721) . offset can most straightforwardly be interpreted as arising from 
while the somewhat gentler removal of a larger-scale reservoir the fading of the SOs relative to their spiral progenitors, with 
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the spread in the relation arising from the diff erent epochs at 
which this fading commenced. Indeed, iBedreeal et all J200d) 
were able to uncover some evidence that the magnitude of each 
SO's offset from the spiral Tully-Fisher relation depends on 
the time since star formation was cut off, as measured by the 
age of its stellar population, just as this picture would predict. 
However, the evidence is still uncomfortably circumstantial, as 
we have had to assume that the ancient progenitors of these 
SOs respect the same Tully-Fisher relation that we see in the 
nearby spirals, and that the difference between the relations 
for SOs and spirals arises from luminosity evolution: the off- 
set could also arise because these systems have fundamentally 
different mass properties, shifting the relation in rotation speed. 
Even with these assumptions, we can only tie spirals and SOs 
together on a statistical basis, as we do not have a direct con- 
nection between any SO and the properties of its individual pro- 
genitor. 

Ideally, we would like to find some historical record pre- 
served in the current properties of an SO that tells us about its 
luminosity during its early life as a spiral, and that will not 
have been defaced by the transformation process or any sub- 
sequent evolution of the galaxy. Fortunately, just such a record 
exists in the globular cluster (GC) population of the SO. As 
we will see below, GCs have a reasonably well defined spe- 
cific frequency (number per unit galaxy V-band luminosity) 
for spiral galaxies, so the number of GCs provides a useful 
proxy for a spiral galaxy's luminosity. Further, the transforma- 
tion from spiral to SO is unlikely to alter the number of GCs 
significantly: the hydrodynamic processes that strip out the gas 
from the spiral will not have any impact on these collisionless 
stellar clusters, and the conversion is likely to be sufficiently 
benign that we are unlikely to lose clusters or produce any 
new ones as we find in m ore dramatic phenomena like mergers 
lAshman & Zepfi n"998). In addition, the GCs are sufficiently 
old that their own passive fading will be very slow, so the ob- 
served number counts should not decrease significantly during 
the transformation. Thus, the number of GCs in an SO offers a 
reasonably robust indicator of its spiral progenitor's luminos- 
ity, so the specific frequency of GCs in such a system provides 
a direct measure of the ratio of its progenitor's luminosity to 
its current luminosity, which we can compare to the less direct 
indicators of fading such as those provided by the Tully-Fisher 
relation. 

In this paper, we carry out such an analysis, comparing 
the specific frequencies of GCs in SOs to those in spirals in 
Section|2]to see how much these systems must have faded, and 
going on to look for evidence that different amounts of fad- 
ing might be related to different epochs of transformation in 
Section[3] Conclusions are presented in Section|4] 

2. The specific frequency of globular clusters in 
SOs and spirals 

The specific frequency of GCs in a galaxy, S n, is usually de- 
fined relative to the galaxy's V-band luminosity, normalized to 
My = -15: 



Table 1. Spiral and SO absolute magnitudes and global globular 
cluster specific frequencies from Ash man & Zepfl dl998l) 
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where N ( is the total number of GCs and M v is 
the total absolute ma gnitude of the underlying galaxy 
jHarris & van den Berghl ll981). Since S n x N t /Ly, as dis- 
cussed above, we would expect S n to grow as spiral galaxies 
fade into SOs, if, indeed, they do. Thus, by comparing the spe- 
cific frequency of GCs in spirals and SOs, it should be possible 
to check whether S n is larger in SOs than in spirals as this sce- 
nario requires, and to estimate by how much on average spirals 
fade during this transformation, for comparison with less direct 
indicators like the Tully-Fisher relation. 

Unfortunately, data sets of GC specific frequencies in 
galaxies with different morphologies do not exist for any 
reasonably-complete sample of galaxies. However, two entirely 
independent compilations have been constructed, so we can 
make some assessment of the likely systematic uncertainties 
that may arise from any selection biases by comparing the re- 
sults from these two data sets. 

First, we have analyzed the values published by 
lAshman & Zepfl ( fl998l) . excluding dwarf galaxies {My > -17) 
and cDs. The relevant data for spirals and SOs are summarized 
in Tabled In Figure^a), we plot the cumulative distributions 
of GC specific frequencies determined for spirals and SOs. For 
completeness, we also show the distribution for elliptical galax- 
ies, but these have no direct bearing on the current analysis. 
It is immediately apparent that SOs have, on average, signifi- 
cantly larger specific frequencies than spirals. A Kolmogorov- 
Smirnov test rejects the hypothesis that spirals and SOs have 
the same S # distribution at the 90% confidence level, and we 
find that for this sample 



(s s N p ) 



= 3.1 + 1.1 



(2) 



s N = Mi(r a4(M £ +15) 



(i) 



where (5'^ p ) and (S^) are the unweighted average GC specific 
frequencies for spirals and SOs respectively. The uncertainties 
on the specific frequencies were calculated from the standard 
error in the mean, and these uncertainties were propagated to 
provide the quoted error in equation|5] 
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Fig. 1. (a) Cumulative distribution of global GC specific frequencies for spiral, SO and elliptical galaxies from Ash man & Zepfl 
(1998). (b) Cumulative distribution of local GC specific frequencies for spiral, SO and elliptical galaxies from Kundu & Whitmore 
(2001a,b) and lGoudfrooii et all J2003k 



We have construct ed the second Sn d a ta set using the 
results published bv iKundu & Whitmorel ll2001al) for el- 
lipti cal galaxies, IKundu & Whitmorel (12001b) for SOs, and 
iGoudfrooii et all J2003I) for spiral galaxies. Again, dwarf 
galaxies have been excluded. The adopted data for spirals and 
SOs are presented in Table [2] One complication here is that 
Kundu and Whitmore published mainly "local" S n values, cor- 
responding to one WFPC2 field, while Goudfrooij et al. give 
global values comparable to those described above. However, 
Kundu & Whitmore (2001a,b, hereafter KW) give both local 
and global specific frequencies for 14 galaxies, so it is pos- 
sible to calibrate the difference between these quantities, and 
how they might depend on the galaxies' absolute magnitude. A 
least-squares fit to these data yields a calibration of 



c local 
^ global 



= 5.71 +0.25M, 1 , 



(3) 



with a scatter of only 0.18. Since we are primarily interested 
in population average values of S and we are looking at vari- 
ations between populations that amount to factors of several, 
this ~ 20% uncertainty is completely negligible. Ideally, we 
would use this transformation to convert the local values from 
KW into global values. However, the calibration of this relation 
depends almost entirely on galaxies with absolute magnitudes 
Mv < -20, and a number of galaxies in the KW sample are 
fainter than this limit, so such a transformation would have sig- 
nificant uncertainties. F ortunately, all the spiral galaxies in the 
IGoudfrooii etail J2003I) sample lie at brighter magnitudes, so 
we can use the relation to transform their quoted global values 
into local values, as shown in Table [2J since we are only in- 
terested in the ratio of specific frequencies, it does not matter 
whether the local or global values are employed as long as we 
use the same in both samples. In figure^b) we plot the cumu- 
lative distributions of GC local specific frequencies for spirals, 
SOs and ellipticals from this second compilation. Again, SOs 
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Table 2. Absolute ma gnitudes and globular cluster specific frequencies from iKundu & Whitmorel (12001 bl) for SOs, and 
iGoudfrooii et alJj200 3) for spiral galaxies, with global colours for SOs where available. 

Name My 5f M 5^ cal (U - B) (B - V) (V - R) (V - 1) 



Spiral Galaxies 



NGC 3628 


-21.0 


1.9 = 


t0.2 


1.0 ±0.1 










NGC 4013 


-20.8 


1.1 = 


t0.3 


0.6 + 0.2 










NGC 4517 


-21.6 


0.6 = 


t0.2 


0.2 ±0.1 










NGC 4565 


-21.5 


0.6 = 


t0.2 


0.2 ±0.1 










NGC 4594 


-22.2 


1.7 = 


t0.6 


0.4 ±0.2 










IC5176 


-21.1 


0.5 = 


tO.l 


0.2 ±0.1 










NGC 7814 


-20.5 


0.7 = 


t0.2 


0.5 ±0.2 










Average 








0.4 ±0.1 











SO Galaxies 



NGC 524 


-22.4 


1.1 ±0.4 


- 


- 


- 




- 


NGC 2768 


-22.0 


1.2 ±0.4 


- 




- 




- 


NGC 6861 


-21.8 


3.6 ±1.6 


0.64 ± 0.03 


1.04 • 0.03 


0.61 ±0.03 


1 


.31 ±0.03 


NGC 6703 


-21.8 


1.2 ±0.8 


- 


- 


- 




- 


NGC 1553 


-21.5 


0.5 ±0.1 


0.49 ± 0.03 


0.95 ± 0.03 


0.57 ± 0.03 


1 


.19 ±0.03 


NGC 474 


-21.4 


0.7 ± 0.5 


- 




- 




- 


NGC 3115 


-21.3 


1.3 ±0.1 


0.59 ± 0.03 


0.99 ± 0.03 


0.61 ±0.03 


1 


.27 ± 0.03 


NGC 1332 


-21.2 


2.2 ± 0.7 


0.60 ± 0.03 


1.03 ±0.03 


0.60 ± 0.03 


1 


.30 ± 0.03 


NGC 3414 


-21.0 


1.6 ±0.6 


0.55 ± 0.01 


0.87 ± 0.01 


0.70 ± 0.01 




- 


NGC 4459 


-20.9 


0.9 ± 0.3 


0.47 ± 0.01 


0.96 ±0.01 


0.55 ± 0.01 


1 


.22 ± 0.02 


NGC 1201 


-20.8 


1.1 ±0.5 


0.54 ± 0.03 


0.97 ± 0.03 


0.57 ± 0.03 


1 


.22 ± 0.03 


NGC 1400 


-20.6 


2.9 ±1.1 


0.61 ± 0.03 


1.03 ±0.03 


0.65 ± 0.03 


1 


.37 ± 0.03 


NGC 3607 


-20.4 


1.3 ±0.4 


- 


- 


- 




- 


NGC 4203 


-20.2 


1.5 ±0.5 


- 


- 


- 




- 


NGC 2902 


-20.2 


0.3 ± 0.4 


— 


— 


— 




— 


NGC 3489 


-19.6 


1.4 ±0.7 












NGC 4379 


-19.6 


0.6 ±0.4 


0.39 ± 0.01 


0.85 ±0.01 


0.61 ±0.01 


1 


.24 ± 0.02 


NGC 1389 


-19.5 


0.5 ±0.4 


0.43 ± 0.01 


0.92 ±0.01 


0.56 ± 0.01 


1 


.16 ±0.02 


NGC 3056 


-18.9 


0.6 ±0.7 












NGC 3156 


-18.9 


0.4 ± 0.4 












IC 3131 


-18.9 


0.2 ±0.4 












NGC 1375 


-18.8 


0.5 ± 0.7 


0.31 ±0.03 


0.81 ±0.03 


0.51 ±0.03 


1 


.06 ± 0.03 


VCC 165 


-18.5 


0.0 ±2.1 












NGC 3599 


-18.4 


1.2 ±0.9 












NGC 2328 


-18.4 


0.0 ±1.4 












NGC 4431 


-18.4 


0.5 ± 0.7 












IC 1919 


-18.2 


1.1 ±1.7 












NGC 1581 


-18.2 


0.2 ± 0.9 


0.20 ± 0.03 


0.73 ± 0.03 


0.46 ± 0.03 


1 


.06 ± 0.03 


ESO 358-G059 


-18.0 


0.6 ±1.3 












NGC 3870 


-17.9 


0.8 ±1.3 












NGC 3115-DW1 


-17.9 


6.8 ± 2.4 












ESO 118-G034 


-17.6 


0.3 ± 2.0 












IC 3540 


-17.3 


3.2 ±4.2 












NGC 4150 


-16.6 


1.7 ±1.8 












Average 




1.2 ±0.2 













have, on average, significantly larger specific frequencies than 
spirals (once again, a Kolmogorov-Smirnov test rejects the hy- 
pothesis that spirals and SOs have the same S n distribution at 
the 90% confidence level, giving a rejection from the combined 
data sets that is significant at the 99% level). For this second 
dataset, we find that the ratio of local specific frequencies is 



Given that equations|2]and|4]were derived from completely dif- 
ferent data sets derived using different instruments using differ- 
ent analysis techniques applied by different authors, the agree- 
ment between the two values suggests that this result is reason- 
ably robust, notwithstanding the incompleteness of the sam- 
ples. It would appear that spirals have a specific frequency of 
GCs that is a factor of ~ 3 lower than that in SOs, with the fairly 
simple interpretation that the latter have held onto their GCs but 
faded by a factor of ~ 3 in transforming from the former. 
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Is this fa ctor plausible? We l l, pop ulation synthesis mod- 
els based on iBruzual & Charloil d2003l) indicate that a galaxy 
with an initially steady star formation rate will fade by a fac- 
tor of ~ 3 in the V band in ~ 3 Gyr after complete truncation 
of its star formation. If the star formation decays more slowly 
(as through the strangulation mechanism), such fading could 
take ~ 4-6 Gyr Fsee lShiova et alJd2004l) for some model exam- 
ples]. These timescales are comparable t o the look-back time 
to z ~ 0.5, which is where Dressi er et al.l (11997) found signs of 
the transition from spiral-dominated to SO-dominated clusters. 
However, we would not expect all galaxies to undergo the tran- 
sition at the same time, and presumably some fraction of the 
significant scatter in Sn for SOs can be attributed to the differ- 
ent amount of fading that has taken place for different systems 
that underwent this transformation at different times, a point to 
which we return in Section[3] Nonetheless, it is heartening that 
the difference in the mean value of S n seems entirely consistent 
with the characteristic look-back time at which the transforma- 
tion is believed to have predominantly occurred. 

The difference also matches up to the results from the 
analysis of the Tuliy -Fisher relation. As mentioned above, 
iBedregal et al.l d200 6) found a ~ 1.4 magnitude offset between 
the B-band Tully-Fisher relation for SOs and that for spirals. 
Using the same passively-evolving stellar population models as 
described above, we find that such a fading in the Z?-band cor- 
responds to a fading in the V-band of ~ 1 .2 magnitudes, or a 
factor of ~ 3 decrease in V-band luminosity, just as found in the 
analysis of globular clusters. Although most reassuri ng, the ex- 
act ag reement must be to an extent coincidentak Bedregal et alJ 
(2006) found significant uncertainty in the offset in the Tully- 
Fisher relations depending on the exact calibration adopted, so 
the factor of ~ 3 predicted by the Tully-Fisher analysis could 
easily be a factor of two or four. The lack of a complete sam- 
ple of SOs upon which to carry out this test also means that 
some residual uncertainty must remain in the absolute value of 
the mean fading of the population, due to unknown systematic 
effects in the way that the sample was selected. 

Thus, although the mean degree of fading as derived from 
the globular cluster specific frequency is entirely consistent 
with the values predicted by population synthesis and those 
derived independently from the Tully-Fisher relation, the case 
would be even stronger if there were a more subtle test that goes 
beyond simple absolute shifts in the mean of the population. 

3. Correlation between globular cluster specify 
frequency and galaxy colours 

Happily, there is a differential test that we can apply to the GC 
specific frequency data that overcomes these uncertainties in 
both absolute magnitude and sample selection. As an individual 
system fades from spiral to SO, the specific frequency of GCs 
will go up as we have discussed above, and the colours of the 
galaxy should redden as the stellar population ages. Thus, we 
should expect to see a correlation between the measured values 
of S n and the SOs' colours, with different values corresponding 
to the degree to which the SO has evolved away from its spiral 
progenitor. Further, we can compare the form of this relation 



to that predicted by the population synthesis calculations to see 
whether the correlation matches the model predictions in detail. 

Unfortunately, the necessary data do not exist for all of 
the SO galaxies for which Sn values are available. The evo- 
lution in colour in such a fading population is relatively mod- 
est, so accurate photometry is required, and we are interested 
in the galaxies' global stellar populations, so need spatially- 
resolved photometry rather than just central values for colours. 
Such high-quality extended data are not rea dily available for 
many galaxies, but using the photometry of IPoulain & Nietol 
(1994) a nd the highes t -quali ty data from the compilation of 
IPrugniel & Heraudeaul dl998l) . we have been able to bring to- 
gether reliable m ulti-colour photometry for 12 of the non-dwarf 
SO galaxies from Kundu & Whitmore (2001b), which we sum- 
marize in Table 13 These 12 galaxies appear to be a fair sub- 
sample of the original sample, in that there is no obvious bias 
in magnitude or the values of S n- To obtain a measure of the 
global populations in these galaxies without unduly compro- 
mising the photometric accuracy, we have calculated the inte- 
grated colours for these systems within an aperture that corre- 
sponds t o the ij-band half-light ra dius of each galaxy, as tabu- 
lated in lde Vaucouleurs et alJ d!976l) . 

Figure |2] shows the derived global U - B, B - V, V - R 
and V — I colours for these 12 SO galaxies plotted against their 
GC specific frequencies. A clear correlation is observed in the 
sense predicted by the fading spiral hypothesis: the larger the 
GC specific frequency in an SO, the redder its stellar popula- 
tion. 

Further, the relation is of exactly the right amplitude. 
Figure |2 shows the predictions of population synthesis codes 
for a fading model stellar population with solar metallicity and 
a standard initial mass function, which ceased forming stars 
more than ~ lGyr ago - after this time, the colour evolu- 
tion is not very sensitive to the detailed history of the star 
formation that took place before star formation ended, so the 
model predictions are fairly robust. To assess the uncertainty 
in the predict ions, we compu ted models using the codes of 
both lBruzual & Chariot! d2003l) and lworfhevl dl994l) . but, as is 
clear from Figure |2] these uncertainties are also fairly negli- 
gible compared to the error bars on the data. In these calcula- 
tions, the variation in colour with age comes directly from the 
models, and the changes in S n come from the corresponding 
variations in the model values of the V-band magnitude, using 
equation and following the basic precept of this paper that 
the number of GCs does not change with time. The jc-axis zero- 
point of the models is fixed so that for t he typical V - I colou r 
of a star-forming spiral, (V - I) ^ 1.1 (Fukugi ta et alJll995l) . 
the models predict S n — 0.4, which is the average value of 
S j^ cal for the star-forming spirals in the lGoudfrooii et ail (|2003 ) 
sample (see Table |2j. Once we have made this calibration to 
the local spirals, the zero-points of the jt-axes in all the other 
plots, the amplitude of all the colours' variations with S n, their 
slopes, and even their curvatures are entirely fixed by the stel- 
lar evolution models, and, as Figure |2] attests, the match of all 
these properties to the observations is remarkably good. 
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Fig. 2 . Global U - B, B - V, V — R and V - I colours versus globular cluster specific frequency for SO galaxies from 
Kun du & Whitmord feOOlbl). The solid and and dashed lines correspond to population synthesis model predictions of fading 
stellar populations from lBruzual & Char iot (2003) and Worthey ( 1994) respectively. 



4. Conclusions 

The classic problem in trying to address questions related to 
galaxy evolution is that our snapshot view of the Universe 
means that we never get to view the process directly. Thus, 
although it seems quite plausible that spiral galaxies fade to 
become SOs when their interstellar media are stripped out and 
star formation ceases, there is little direct evidence to this ef- 
fect. Here, we finesse this issue by making use of the histori- 
cal record imprinted in a galaxy via its GC population, which 
should still be "readable" even after the transformation pro- 
cess. In particular, if the number of globular clusters is taken 
as a measure of the progenitor spiral galaxy's luminosity, and 
this number does not change as the galaxy passively evolves, 
then we have been able to show that SOs must have faded by an 
average factor of approximately three during this process. This 
result matches very neatly with the values predicted by stellar 
population synthesis models and those found by comparing the 
Tully-Fisher relations of spirals and SOs. 



Further, using the amount by which the stellar population 
has reddened as a measure of the time since the transformation 
into an SO began, we have been able to show that individual 
galaxies are at different stages along this evolutionary track: the 
colours of these galaxies correlate with the specific frequency 
of GCs in exactly the way that one would expect if their star 
formation shut down at different times in the past. Thus, we are 
at the point of being able to follow the histories of individual 
SO galaxies with the detail necessary to say which underwent 
their life-changing transformations first. 

The next logical step in writing the life histories of these 
galaxies is to obtain uniform samples of spectra of the quality 
nece ssary to derive stellar popul ation ages for these systems 
(e.g. iKuntschner & Davieslll998i) . Not only would these data 
resolve any remaining ambiguity by lifting the degeneracy be- 
tween age and metallicity effects in the interpretation of the 
broad-band colours used here, but they would also allow us to 
put much more accurate dates to the "birthdays" of SO galaxies. 
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